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ABSTRACT 


Six  events  from  the  Baikal  rift  zone  occurring  from  1971  to  1975  were  ex- 
amined in  a seismic  discrimination  context.  Seismograms  from  ALPA,  LASA,  NOR- 
SAR,  the  HGLP,  the  WWSSN,  and  the  SRO  stations  were  studied  for  source  mechan- 
ism, Ms-m^,  corner  frequency,  pP,  complexity,  and  spectral  ratio.  All  the 
Baikal  events  can  be  identified  as  earthquakes  except  the  November  5,  1976, 
event  which  shows  explosion  characteristics.  The  characteristics  of  the  Baikal 
events  were  compared  to  six  East  Kazakh  explosions  and  the  Aleutian  explosion 
MILROW.  No  evidence  is  found,  from  the  signals  studied  here,  for  a zone  of  high 
attenuation  under  the  Baikal  rift.  The  discriminants  found  to  be  useful  in 
this  study,  with  the  possible  exception  of  short-period  spectral  ratio  which 
is  rather  unreliable,  would  not  give  correct  discrimination  for  a suitable 
shot  array.  However  there  is  a full  magnitude  range  of  separation  in  M^m^ 
in  this  region,  suggesting  that  a suitable  shot  array  might  be  difficult  to 
devise. 
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INTRODUCTION 


Discrimination  parameters  have  been  applied  in  detail  at  only  a few  nu- 
clear test  sites  and  earthquake  regions  and  the  wide  variability  of  path-re- 
lated factors  ai.ecting  discrimination  is  not  yet  fully  revealed.  This  study 
is  part  of  a series  which  will  extend  these  discrimination  studies  to  other 
regions  of  shallow  earthquakes.  This  report  examines  six  such  events  in  the 
Baikal  rift  zone,  which  is  shown  in  Figure  1.  Most  of  the  earthquakes  in  the 
Baikal  rift  zone  are  probably  associated  with  normal  faults  bordering  graben 
structures  in  the  rift  zone.  Actual  seismogram  analysis  revealed  that  the 
events  chosen  were  too  small  to  determine  their  exact  fault  plane  solutions. 

We  determine  average  Mg  and  m^  for  the  events  in  this  study  and  apply  other 
common  discrimination  parameters  such  as  first  motion,  comer  frequency,  com- 
plexity, and  short-period  spectral  ratios  to  the  data.  We  examine  source  and 
propagation  affects  on  the  principle  discrimination  phases,  P and  LR.  Results 
for  the  Baikal  events  will  be  compared  to  six  East  Kazakh  explosions  and  the 
Aleutian  explosion  MILROW. 


* NEIS  avants  (1971-1975)  chasan  far  this  study 

♦ additional  NEIS  apicantars  (1971  1975) 

^ saismic  statian  locatians 

mr  normal  faults 
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Figure  1. 


MONGOLIA 


Map  of  the  Baikal  rift  zone  showing  faults  observed  on  satellite 
photographs  and  NEIS  epicenters  from  January  1971  through  Febru- 
ary 1975. 


TECTONIC  SETTING 


General  Features 

Figure  1 is  a map  of  the  Baikal  rift  zone,  a system  of  grabens  extending 
from  Lake  Baikal  to  about  2000  km  northeast.  The  grabens  are  a geologically 
young  feature  situated  in  a Paleozoic  fold  belt.  The  Baikal  rift  zone  is  a 
region  of  normal  faulting  and  crustal  extension  which  may  be  a remote  expres- 
sion of  the  stress  system  resulting  from  the  Indian-Eurasian  collision  (Mol- 
nar  and  Tapponnier,  1975).  According  to  U.S.S.R.  seismic  studies  (Solonenko, 
1968),  seismicity  is  confined  to  the  crust. 

Positive  isostatic  gravity  anomalies  follow  the  rift  zone.  The  best  ex- 
planation for  these  anomalies  is  a raised  Moho  boundary  under  the  rift  zone 
which  overlies  mantle  of  anomalously  low  density  (Zorin,  1971).  Heat  flow 
measurements  in  Lake  Baikal  are  higher  than  normal  (Lubimova,  1969) , confirm- 
ing the  magneto-telluric  detection  of  a highly  conducting  layer  in  the  upper 
mantle.  Abnormal  heating  of  the  upper  mantle  beneath  the  rift  zone  may  be 
causing  the  tectonic  activity  in  this  region. 

Source  Mechanisms  of  Earthquakes 

Most  of  the  earthquakes  in  the  rift  zone  are  associated  with  the  normal 
faults  which  comprise  the  boundaries  of  the  grabens  (Solonenko,  1968;  Molnar 
et  al.,  1973). 


Molnar,  P.,  and  P.  Tapponier  (1975).  Cenozoic  tectonics  of  Asia:  effects  of 
a continental  collision.  Science.  189,  419. 

Solonenko,  V.  P.  (1968).  Regional  seismic  zoning  of  the  U.S.S.R.  - Eastern 
Siberia,  in  Seismic  Zoning  of  the  U.S.S.R.,  edited  by  S.  V.  Medvedev, 
English  translation,  Keter  Publishing  House,  Jerusalem. 

Zorin,  Y.  (1971).  Recent  Structure  and  Isostasy  of  the  Baikal  Rift  Zone  and 
Surrounding  Territory.  Nauka,  Moscow. 

Lubimova,  E.  (1969).  Heat  flow  patterns  in  Baikal  and  other  rift  zones,  Tec- 
tonophysics.  8,  457. 

Molnar,  P.,  T.  Fitch,  and  F.  Wu  (1973).  Fault  plane  solutions  of  shallow 
earthquakes  and  contemporary  tectonics  in  Asia,  Earth  and  Planetary 
Science  Letters,  19,  101. 
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NEIS  epicenters  for  events  in  this  area  from  January  1971  through  February 
1975  are  shown  in  Figure  1 overlaying  faults  which  were  inferred  from  our  sat- 
ellite photomosaic  of  the  area.  Event  b occurred  in  1976  and  was  u special 
addition  to  the  data  set.  It  can  be  seen  that  the  earthquakes  are  near 
the  fault  traces  except  for  event  6.  Studies  of  earthquake  focal  mechanisms 
in  the  Baikal  region  have  shown  that  the  tensile  stresses  are  horizontally  ori- 
ented, perpendicular  to  the  trend  of  the  grabens  (Vvedenskaya,  1961). 

a 

Velocity.  Density,  and  Attenuation 

Seismic  studies  in  the  Baikal  region  report  an  average  crustal  thickness 
of  36  to  38  km  in  the  rift  zone  and  36  to  46  km  in  the  surrounding  region 
(Beliayevsky  et  al.,  1968).  Geophysical  studies  indicate  that  the  upper  man- 
tle under  the  rift  zone  has  low  velocities,  low  density,  and  possibly  high 
attenuation  (Kodean,  197b).  Areas  of  high  attenuation  have  also  been  observed 
southwest  of  the  rift  zone  (Vinnik  and  Godzikovskaya,  1972,  1973).  Table  I 
shows  a typical  velocity  profile  in  the  rift  zone  derived  from  deep  seismic 
sounding  studies  (Puzvrev  et  al.,  1973). 

The  Baikal  region  has  been  compared  with  the  African  rift  zone  and  has 
many  features  similar  to  the  Basin-Range  Province  in  the  western  United  States, 
also  a Cenozoic  extensionnl  region.  Of  special  interest  in  this  study  is  the 
degree  of  attenuation  for  P waves  under  the  Baikal  source  region.  High  atten- 
uation has  been  linked  to  positive  travel-time  residuals  in  the  western  United 


Vvedenskaya,  A.  (1961).  Special  features  of  the  stressed  state  in  foci  of 
earthquakes  in  the  Baikal  region,  Izvestla,  Geophys.  Ser..  432. 

Beliayevsky,  N.,  A.  Borisov,  1.  Valvovsky,  and  Y.  Schukin  (1968).  Transcon- 
tinental crustal  sections  of  the  U.S.S.R.  and  adjacent  areas,  Canadian 
Journal  of  Earth  Sciences,  5,  1067. 

Rodean,  H.  (1976).  Seismic  yield  verification  and  a regional  vs.  m^  anomaly, 

Lawrence  Livermore  Laboratory,  UCID-17006,  16  pp. 

Vinnik,  L.,  and  A.  Godzikovskaya  (1972).  Sounding  of  the  Earth's  mantle  by  the 
method  of  seismically  conjugate  points,  Izvestla,  Earth  Physics.  No.  10, 

656. 

Vinnik,  L.,  and  A.  Godzikovskaya  (1973).  Lateral  variations  of  the  absorption 
by  the  upper  mantle  beneath  Asia,  Izvestla.  Earth  Physics.  No.  1,  1. 

Puzyrev,  N. , M.  Mandelbaum,  S.  Krylov.  B.  Mishenkin,  G.  Krupskaya,  and  G. 

Petrlck  (1973),  Deep  seismic  investigations  in  the  Baikal  rift  zone, 

Tectonophysics.  20.  85. 
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TABLK  I 


States  (Der,  1976);  and  In  consideration  of  the  similarities  between  the  Baikal 
rift  zone  and  the  western  United  States,  we  Investigated  residuals  for  seismic 
stations  located  In  this  rift.  Only  four  which  regularly  report  to  the  ISC 
(International  Selsmologlcal  Center)  are  relevant:  ZAK,  BOD,  TUP  and  YAK. 

These  locations  are  shown  in  Figure  1;  note  that  they  are  not  In  the  rift  zone 
proper.  The  1974  ISC  bulletins  were  gleaned  for  travel-time  data  at  these 
stations  for  large  earthquakes.  An  effort  was  made  to  sample  somewhat  evenly 
the  teleselsmic  distance  range  and  the  full  range  of  back  azimuths.  The  re- 
siduals (relative  to  Jef f reys-Bullen  tables)  are  plotted  in  Figure  2.  In  gen- 
eral, the  residuals  tend  to  negative  times,  implying  higher  upper  mantle  or 
crustal  velocities  for  media  under  the  stations.  Even  subsets  for  each  sta- 
tion along  that  range  of  azimuth  which  would  enclose  rays  emerging  up  under 
the  Baikal  rift  zone  do  not  apparently  show  a positive  residual,  although  the 
data  is  probably  too  sparse  and  the  stations  are  not  well  enough  situated  to 
draw  definite  conclusions. 


Der,  Z.  A.  (1976).  On  the  existence,  magnitude  and  causes  of  broad  regional 
variations  in  body-wave  amplitudes,  SDAC-TR-76-8,  Teledyne  Geotech,  Alex- 
andria, Virginia. 

I 
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Figure  2.  Selected  ISC  residuals  for  stations  in  or  near  Baikal  rift  zone. 


Selected  ISC  residuals  for  stations  in  or  near  Baikal  rift 
zone. 


Selected  ISC  residuals  for  stations  in  or  near  Baikal  rift 


Figure  2 (cont.)  Selected  ISC  residuals  for  stations  in  or  near  Baikal  rift 


DATA 


Earthquake  Selection 

For  this  study  shallow-focus  earthquakes  were  selected  from  the  Baikal 
rift  zone.  The  NEIS  list  locates  the  large  majority  of  the  earthquakes  in  the 
Baikal  rift  zone  within  the  crust.  The  earthquakes  chosen  were  limited  to  the 
years  1971  to  1975  so  that  we  could  utilize  the  data  from  the  large  seismic 
arrays  and  the  HGLP  network.  We  selected  a total  of  six  events  as  listed  in 
Table  11  and  shown  in  Figure  1.  Again,  we  note  that  event  6 was  in  197b  and 
added  late  to  the  data  set.  These  particular  six  events  were  chosen  on  the 
basis  of  relatively  low  M^,  compared  to  m^,  as  measured  by  Rayleigh  waves 
recorded  at  NORSAR.  Table  II  also  lists  the  depths  determined  from  our  pP 
observations  at  LASA,  NORSAR,  and  the  average  of  selected  WWSSN  stations. 

All  of  the  selected  events  are  located  within  the  crust. 

Explosion  Selection 

The  Baikal  events  were  compared  to  explosions  in  East  Kazakh,  USSR,  which 
lies  to  the  southwest  of  the  Baikal  rift  zone  and  to  the  Aleutian  explosion 
MILROW.  We  selected  a total  of  seven  explosions  as  listed  in  Table  III  for 
this  comparison. 

Seismic  Stations 

Digital  data  from  the  three  arrays  ALPA,  LASA,  and  NORSAR  and  from  the 
available  HGLP  stations  and  film  data  from  selected  WWSSN  stations  were  assem- 
bled for  these  events.  The  available  LASA  and  NORSAR  short-period  P recordings 
are  shown  in  Figure  3.  We  have  used  the  AO  subarray  from  LASA  and  the  C3  sub- 
array from  NORSAR  as  in  our  previous  reports  of  this  series  in  order  to  avoid 
signal  loss  at  high  frequencies  due  to  full-array  beaming.  WWSSN  stations 
were  selected  on  the  basis  of  magnification  and  proximity  to  the  region  of 
study,  and  it  is  unlikely  that  stations  not  used  here  would  significantly  add 
to  the  data  base  for  these  Baikal  events.  Since  WWSSN  data  was  not  available 
at  the  time  of  this  study  for  event  6,  we  examined  SRO  data  for  this  event  in- 
stead of  WWSSN  data. 
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TABLE  III 

NEIS  Parameters  for  East  Kazakh  and  MILROW  Explosions 


EPICENTER 


Date 

Origin  Time 

LAT  (aN) 

LONG  (°E) 

JEb 

Region 

2 Oct  1969 

22  06  00.0 

51.418 

179.181 

6.5 

Aleutians 

22  Mar  1971 

04  32  57.8 

49.744 

78.185 

5.8 

East  Kazakh 

29  Nov  1971 

06  02  57.1 

49.758 

78.126 

5.5 

East  Kazakh 

30  Dec  1971 

06  20  57.7 

49.750 

78.130 

5.8 

East  Kazakh 

10  Feb  1972 

05  02  57.3 

49.986 

78.886 

5.5 

East  Kazakh 

25  Jun  1974 

03  56  57.6 

49.889 

78.115 

4.8 

East  Kazakh 

07  Dec  1974 

05  59  56.9 

49.908 

77.648 

4.8 

East  Kazakh 

SIGNAL  ANALYSIS 


Source  Effects 

We  wanted  initially  to  identify  a source  mechanism  for  each  event  so  that 
we  could  predict  the  radiation  pattern  for  body-wave  and  surface-wave  phases. 
Determination  of  the  fault  planes  for  earthquakes  whose  magnitudes  are  below 
6.0  is  generally  not  reliable  with  teieseismic  data,  and  all  our  earthquakes 
have  ra^  less  than  6.  First  motions  were  determined  when  possible  from  the 
data  we  have  collected,  and  ISC  bulletin  data  was  available  for  events  1 through 
5.  We  have  plotted  all  the  first  motions  for  the  Baikal  events  in  Figure  4.  No 
clear  first  motions  were  available  for  event  6.  Unfortunately,  the  short-period 
data  were  not  consistent,  and  it  is  impossible  on  the  basis  of  the  plots  alone 
to  determine  individual  fault  planes  because  dilatational  and  compressional 
first  motions  do  not  separate  well.  However,  taken  together,  there  is  a domi- 
nance of  dilatational  first  motion  in  the  central  part  of  the  focal  sphere 
plots;  this  would  result  from  normal  faulting  and  is  in  agreement  with  the 
body  of  known  focal  mechanisms  for  this  region. 

Corner  frequencies  and  seismic  moments  for  our  Baikal  events  have  been 
estimated  from  short-period  LASA  and  NORSAR  spectra.  The  spectra  are  shown  in 
Figure  5.  They  are  from  the  phased  beams  of  the  AO  subarray  at  LASA  and  C3 
subarray  at  NORSAR.  The  sample  length  was  6.4  seconds.  The  signals  have  been 
tapered,  and  the  instrument  response  was  removed  from  the  signal  and  noise  spec- 
tra. Attenuation  was  removed  from  the  spectra  by  multiplying  by  the  factor 

-2  -3 

exp  (rtf  t*)  with  a t*  of  .59  for  a go  source  model  and  .40  for  a go  source  model 

-2  -3 

at  LASA  and  t*  of  .18  for  w and  .00  (no  attenuation)  for  co  source  models 

at  NORSAR.  The  basis  for  these  t*  values  will  be  shown  later  in  this  report. 

Corner  frequencies  were  estimated  with  the  assumption  of  complete  stress 
-2  -3 

drop  and  a w or  u asymptotic  relation  at  high  frequencies.  Many  of  the 
estimates  here  of  corner  frequency  and  long-period  level  are  tenuous,  and 
those  which  were  made  with  little  confidence  are  represented  by  dashed  lines 
in  Figure  5.  Occasionally  we  computed  spectra  for  256-second  windows  for 
those  cases  where  the  signal-to-noise  ratio  was  adequate  at  lower  frequencies, 
but  for  these  cases  the  shape  of  the  spectra  was  sufficiently  similar  to  those 
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Figure  5.  LASA  AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Saikal 
events  with  instrument  response  and  attenuation  removed. 
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LASA  AO  and  NORSAR  C3  subarray  spectra  of  P waves  from 
Baikal  events  with  instrument  response  and  attenuation 
removed. 
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Figure  5 (cont.)  LASA  AO  and  NORSAR  C3  subarray  spectra  of  P waves  from 
Baikal  events  with  instrument  response  and  attenuation 
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LASA  AO  and  NORSAR  C3  subarray  spectra  of  P waves  from 
Baikal  events  with  instrument  response  and  attenuation 
removed. 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 


FREQUENCY  (Hz) 

LASA  AO  and  NORSAR  C3  subarrav  spectra  of  P waves  from 
Baikal  events  with  instrument  response  and  attenuation 
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Figure  5 (cont.)  LASA  AO  and  NORSAK  C3  subarray  spectra  of  P waves  from 
Baikal  events  with  instrument  response  and  attenuation 


for  the  6.4-second  windows  that  have  not  shown  them.  Note  that  for  some 
cases  the  signal-to-noise  ratio  was  high  down  to  low  frequencies. 
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Seismic  moment  was  calculated  from  the  body-wave  spectral  displacement 
levels  estimated  in  Figure  5 at  low  frequencies  using  the  relation 

u 4ttpoi3d|J2  I 
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where  was  assumed  to  be  unity  due  to  lack  of  firm  knowledge  of  the  focal 
mechanism,  values  of  pa3  from  Table  I appropriate  to  the  source  depth  were 
used,  and  the  divergence  factor  D was  applied  (Ben-Menahem  et  al.,  1965).  A 
graph  of  moment  versus  corner  frequency  is  shown  in  Figure  6.  Our  events  fall 
close  to  or  between  the  constant  stress  drop  lines  of  Hanks  and  Thatcher  (1972), 
indicating  that  our  events  are  of  intermediate  stress  drop. 

We  attempted  to  find  the  LR  radiation  pattern  for  each  event  by  plotting 
the  antilogs  of  the  Mg  values  for  each  event  as  listed  in  Table  IV.  The  re- 
sults from  the  event  with  the  largest  Mg,  event  1,  are  shown  in  Figure  7.  The 
maximum  antilog  of  Mg  for  the  event  plots  on  the  outer  circle  and  all  the  other 
antilog  values  for  the  event  are  scaled  linearly  between  the  center  of  the 
plot  and  the  radius  of  the  circle.  The  observed  distribution  of  amplitudes 
does  not  obviously  conform  to  any  quadripole  radiation  pattern,  which  is  prob- 
ably due  to  varying  propagation  effects  coupled  with  the  large  epicentral  dis- 
tances and  the  small  magnitude  of  events  studied.  The  scatter  of  data  points 
for  all  the  events  was  so  poor  that  .o  attempt  was  made  to  find  the  best-fit- 
ting LR  radiation  patterns. 

Propagation  Effects 

We  measured  relative  mantle  attenuation  along  the  paths  from  the  Baikal 
events  to  LASA  and  NORSAR,  the  only  two  sites  for  which  we  had  short-period 
digital  data.  Our  measurement  of  attenuation  is  the  factor  t*,  which  is  the 


Ben-Menahem,  A.,  S.  Smith,  and  T.  Teng  (1965).  A procedure  for  source  studies 
from  spectrums  of  long-period  seismic  body  waves.  Bull.  Seism.  Soc.  Am. , 

55,  203. 

Hanks,  T.,  and  W.  Thatcher  (1972).  A graphical  representation  of  seismic  source 
parameters,  J_.  Geophys.  Res. , 77,  4393. 
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Seismic  moment  versus  corner  frequency  for  Baikal  events  from 
LASA  and  NORSAR  P recordings. 


negative  M values  reflect  a noise  measurement  when  no  signal  is  visible 


EVENT  3 


Station 

A 

% 

M 

s 

AAE 

70.1 

4.58 

ATU 

57.4 

5.53 

BUL 

99.4 

5.00 

HLW 

58.7 

5.23 

4.40 

JER 

54.9 

4.80 

4.72 

ROD 

49.0 

5.46 

5.21 

LAO 

76.2 

5.02 

4.50 

MAT 

26.9 

4.24 

4.03 

MSH 

37.2 

5.06 

4.60 

NAI 

79.5 

4.80 

4.92 

NAO 

48.5 

4.21 

5.31 

NIH 

32.9 

4.88 

4.55 

NIL 

31.0 

4.58 

4.91 

NIIR 

43.1 

4.70 

5.04 

POO 

43.1 

4.97 

QUE 

37.2 

4.75 

SHI 

46.0 

5.37 

5.05 

SHL 

29.8 

5.24 

TAB 

44.0 

4.24 

TRI 

57.3 

4.77 

5.51 

UME 

43.4 

4.68 

5.15 

AVERAGE 

4.92 

4.79 

Event  4 


Station 

A 

% 

M 

8 

ALP 

40.5 

-3.52 

ALQ 

79.4 

4.73 

4.51 

CHG 

41.7 

4.57 

3.89 

COL 

40.7 

5.17 

CTA 

78.8 

4.91 

4.27 

EIL 

65.0 

5.76 

4.60 

JER 

63.1 

4.91 

4.30 

KBL 

42.4 

5.36 

4.34 

DIP 

66.0 

4.37 

KON 

50.0 

4.75 

KON 

50.0 

4.55 

LAO 

68.9 

5.48 

NAO 

51.5 

4.96 

4.55 

NIL 

40.8 

4.75 

NUR 

47.0 

4.64 

OGD 

81.4 

4.84 

QUE 

46.9 

4.59 

4.50 

SHL 

38.2 

4.58 

TUI 

74.1 

4.  76 

4.32 

TLO 

74.1 

4.57 

AVERAGE 

4.95 

4.37 
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travel  time  along  the  ray  path  divided  by  the  average  attenuation  factor  Q. 

We  determine  t*  by  first  normalizing  and  then  summing  the  t*=0  spectra 
shown  in  Figure  5 for  each  of  the  two  paths:  Baikal  to  LASA  and  Baikal  to 
NORSAR.  These  averages  are  shown  in  Figures  8 and  9 and  represent  the  aver- 
age relative  ground  displacement  density  spectra  at  the  two  sites.  If  we 

_2 

assume  that  the  source  spectra  has  an  f fall  off  at  high  frequencies,  then 
F(f ) * S.f"2.e_1Tft* 

for  frequencies  above  the  corner  frequency  where  F(f)  is  the  observed  spectrum 
(corrected  for  seismograph  response),  S is  the  source  spectrum  scale  constant, 
and  f is  frequency  in  Hz.  Then  after  taking  the  natural  log  of  this  equation 
ln[F(f ) ] + 2.1n(f)  = -rrt*f  + ln(S) 

If  we  plot  ln[F(f)]  + 2«ln(f)  versus  frequency  for  each  path,  as  shown  in  Fig- 
ures 10  and  12,  the  slope  of  the  graph  is  -irt*.  We  calculated  the  slope  in  the 
frequency  range  1.0  to  2.5  Hz,  where  the  signal-to-noise  ratio  was  highest, 

and  presumably  beyond  the  corner  frequency  so  that  the  high-frequency  asymp- 

-2 

totic  slope  of  f characterizes  the  source  spectrum.  Figure  5 shows  that 
some  of  the  corner  frequencies  are  above  1 Hz,  so  the  data  may  not  be  well 
fitted  by  the  above  equation;  but  if  a higher  frequency  range  were  used,  we 
would  seriously  contaminate  the  average  spectra  with  background  noise.  The 
effect  of  using  spectral  estimates  at  frequencies  below  the  corner  frequency 
is  to  bias  the  slope  t*  toward  a lower  value  than  the  actual  one.  Thus  t* 
results  here  are  not  as  reliable  as  needed  to  confidently  characterize  the 
attenuation  along  the  two  paths.  Similarly,  the  equation  for  a source  spectra 
with  an  f ^ fall  off  is 

ln[F(f ) ] + 3.1n(f)  - -irt*f  + ln(S) 

If  we  plot  the  left-hand  side  of  this  relation  versus  frequency  for  each  path 
as  shown  in  Figures  11  and  13,  the  slope  of  the  graph  is  again  -vt*;  but  t*  is 

_3 

now  based  on  an  assumed  f source  spectrum. 

_2 

The  t*  values  for  an  f source  model  are  .59  + .40  for  the  Baikal  to 
LASA  path  and  .18  + .07  for  the  Baikal  to  NORSAR  path.  The  t*  values  for  an 

_3 

f source  model  are  .40  + .16  for  the  Baikal  to  LASA  path  and  -.00  + .07  for 
the  Baikal  to  NORSAR  path.  Separate  t*  values  were  calculated  for  event  6 
since  it  lies  a few  degrees  to  the  north  of  the  rift  zone  and  because  it  is 


r 
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presumed  to  be  an  explosion  on  the  basis  of  data  to  be  presented  later.  The 
-2 

t*  value  for  an  f source  model  is  -.01  + .19  for  the  path  to  NORSAR  and  for 
-3 

an  f source  model  the  t*  value  for  the  path  to  NORSAR  is  -.20  + .18.  Since 

the  value  of  t*  cannot  be  less  than  zero,  the  t*  value  for  the  path  to  NORSAR 

-2  -3 

is  assumed  to  be  zero  for  both  the  f and  f source  models.  We  add  that 

-3 

there  is  a predominance  of  theoretical  models  which  imply  a f source  spec- 

-2 

trum  rather  than  an  f one  (Geller,  1976).  There  appears  to  be  a real  dif- 
ference in  attenuation  between  the  two  paths.  The  Baikal  to  NORSAR  paths  have 
less  attenuation  than  to  LASA,  as  expected  since  the  former  path  is  mostly 
through  continental  shield  areas  while  the  latter  path  includes  the  anomalous 

upper  mantle  of  the  western  United  States.  Values  of  t*  for  the  Kazakh  explo- 
-2 

sions  for  an  f source  model  (von  Seggern  and  Blandford,  1972)  were  also  cal- 
culated. These  t*  values  are  .16  + .10  for  the  East  Kazakh  to  LASA  path,  .05 
+ .08  for  the  East  Kazakh  to  NORSAR  path,  and  .23  + .10  for  the  MILROW  to  LASA 
path.  As  expected  from  the  reason  just  stated,  the  t*  values  for  the  paths  < > 
NORSAR  are  lower  than  the  values  for  the  paths  to  LASA. 

von  Seggern  and  Sobel  (1975)  have  presented  a method  of  predicting  20- 
second  Rayleigh-wave  amplitudes  by  geometrical  ray  tracing.  Their  work  shows 
that  lateral  changes  in  20-second  LR  phase  velocity  can  cause  major  multipath- 
ing  and  also  large  teleseismic  amplitude  variations  due  to  accumulated  small- 
scale  refraction  effects.  Figures  14  and  15  illustrate  the  paths  of  rays  which 
result  for  events  in  the  eastern  and  western  parts  of  the  rift  zone.  No  at- 
tempt is  made  to  quantitatively  relate  these  patterns  to  observed  LR  amplitudes 
since  we  cannot  correct  the  predicted  amplitudes  for  source  mechanism,  but  we 
can  qualitatively  compare  the  pattern  of  rays  to  residuals  calculated  from 
the  magnitude  data  in  Table  IV.  (The  rays  which  pass  close  to  the  North  Pole 
were  stopped  close  to  that  point  due  to  ray-tracing  numerical  problems  near 


Geller,  R.  J.  (1976).  Scaling  relations  for  earthquake  source  parameters  and 
magnitude.  Bull.  Seism.  Soc.  Am. , 66,  1501. 

von  Seggern,  D.,  and  R.  Blandford  (1972).  Source  time  functions  and  spectra  for 
underground  nuclear  explosions,  Geophys . , 31,  83. 

von  Seggern,  D.,  and  P.  Sobel  (1975).  Experiments  in  refining  M estimates  for 
seismic  events,  SDAC  Report  No.  TR-75-14,  Teledyne  Geotech,  Alexandria, 
Virginia. 
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presumed  to  be  an  explosion  on  the  basis  of  data  to  be  presented  later.  The 
-2 

t*  value  for  an  f source  model  Is  -.01  + .19  for  the  path  to  NORSAR  and  for 
-3 

an  f source  model  the  t*  value  for  the  path  to  NORSAR  is  -.20  + .18.  Since 

the  value  of  t*  cannot  be  less  than  zero,  the  t*  value  for  the  path  to  NORSAR 
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is  assumed  to  be  zero  for  both  the  f and  f source  models.  We  add  that 

-3 

there  is  a predominance  of  theoretical  models  which  imply  a f source  spec- 

_2 

trum  rather  than  an  f one  (Geller,  1976).  There  appears  to  be  a real  dif- 
ference in  attenuation  between  the  two  paths.  The  Baikal  to  NORSAR  paths  have 
less  attenuation  than  to  LASA,  as  expected  since  the  former  path  is  mostly 
through  continental  shield  areas  while  the  latter  path  includes  the  anomalous 

upper  mantle  of  the  western  United  States.  Values  of  t*  for  the  Kazakh  explo- 
-2 

sions  for  an  f source  model  (von  Seggern  and  Blandford,  1972)  were  also  cal- 
culated. These  t*  values  are  .16  + .10  for  the  East  Kazakh  to  LASA  path,  .05 
+ .08  for  the  East  Kazakh  to  NORSAR  path,  and  .23  + .10  for  the  MILROW  to  LASA 
path.  As  expected  from  the  reason  just  stated,  the  t*  values  for  the  paths  to 
NORSAR  are  lower  than  the  values  for  the  paths  to  LASA. 


von  Seggern  and  Sobel  (1975)  have  presented  a method  of  predicting  20- 
second  Rayleigh-wave  amplitudes  by  geometrical  ray  tracing.  Their  work  shows 
that  lateral  changes  in  20-second  LR  phase  velocity  can  cause  major  multipath- 
ing  and  also  large  teleseismic  amplitude  variations  due  to  accumulated  small- 
scale  refraction  effects.  Figures  14  and  15  illustrate  the  paths  of  rays  which 
result  for  events  in  the  eastern  and  western  parts  of  the  rift  zone.  No  at- 
tempt is  made  to  quantitatively  relate  these  patterns  to  observed  LR  amplitudes 
since  we  cannot  correct  the  predicted  amplitudes  for  source  mechanism,  but  we 
can  qualitatively  compare  the  pattern  of  rays  to  Mg  residuals  calculated  from 
the  magnitude  data  in  Table  IV.  (The  rays  which  pass  close  to  the  North  Pole 
were  stopped  close  to  that  point  due  to  ray-tracing  numerical  problems  near 


Geller,  R.  J.  (1976).  Scaling  relations  for  earthquake  source  parameters  and 
magnitude.  Bull.  Seism.  Soc.  Am. , 66,  1501. 

von  Seggern,  D.,  and  R.  Blandford  (1972).  Source  time  functions  and  spectra  for 
underground  nuclear  explosions,  Geophys . ^J. , 31,  83. 

von  Seggern,  D.,  and  P.  Sobel  (1975).  Experiments  in  refining  Mg  estimates  for 
seismic  events,  SDAC  Report  No.  TR- 75-14,  Teledyne  Geotech,  Alexandria, 
Virginia. 
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the  poles.)  Where  the  rays  cross  we  can  expect  either  constructive  interfer- 
ence (large  positive  residuals  as  at  NIL,  OGD,  and  TOL/TLO)  or  destructive 


interference  (large  negative  residuals  as  at  JEK  and  LASA) . The  M residuals 

s 

at  TAB  in  Iran  are  negative,  as  can  be  predicted  from  the  diverging  raypaths 


near  that  site.  The  Mg  residuals  at  KON  and  NORSAR  are  positive  as  can  be 
predicted  from  the  converging  ray  paths  near  those  sites,  and  the  Mg  residuals 
at  QUE  are  about  zero  as  can  be  predicted  from  the  smooth  ray  pattern  near 
that  station.  These  results  are  only  rough  approximations  since  the  velocity 
grid  was  not  well  known  in  some  places  and  attenuation  has  not  been  taken  into 
account.  Also  there  are  many  recorded  amplitudes  which  do  not  correlate  well 
with  the  ray-tracing  predictions.  Stations  which  record  a stable  amplitude 
are  desirable  for  calibrating  a source  region  for  amplitude.  For  events  from 
the  Baikal  region,  such  stations  could  be  located  in  Eastern  Europe,  Western 
USSR,  and  Eastern  China.  The  complicated  LR  ray-tracing  figures  partially  ex- 
plain why  we  were  unable  to  determine  LR  radiation  patterns  from  our  amplitude 


DISCRIMINATION  ASPECTS 
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In  this  and  the  following  sections  we  apply  several  common  discriminants 
between  earthquakes  and  explosions  to  the  Baikal  events.  East  Kazakh  explo- 
sions, and  MILROW.  We  independently  determined  average  Mg  and  m^  values  for 
the  Baikal  events  in  Table  11  using  WWSSN  stations,  HGLP  stations,  the  ALPA, 
LASA,  and  NORSAR  arrays,  and  SRO  data.  The  data  used  in  the  average  estimates 
is  listed  in  Table  IV.  Negative  values  in  the  table  indicate  noise  measure- 
ments which  were  used  as  upper  bounds  for  the  signal  amplitudes.  Magnitudes 
were  computed  according  to  the  formulas 
- log  (A/T)  + B(A) 

Mg  - log  (A/T)  + 1.66  logA  +0.3 

where  A - one-half  the  peak-to-peak  maximum  recorded  amplitude  reduced  to  mvi 
ground  displacements, 

T » period  in  seconds  (restricted  to  17  to  23  sec  for  M calculations), 

A - epicentral  distance  in  degrees, 

B(A)  - Gufenberg-Richter  correction  term  for  surface-focus  P waves. 

We  also  independently  determined  Mg  for  the  East  Kazakh  explosions  in  Table 
III  using  the  NORSAR  array  and  HGLP  stations.  The  data  used  in  the  average 
estimates  for  explosion  Mg  are  listed  in  Table  V,  and  explosion  m^  values  are 
from  the  NEIS  list.  Note  that  the  station  distances  for  both  the  earthquakes 
and  explosions  (Table  IV  and  V)  are  all  teleseismic,  and  therefore  regional- 
ized magnitude  formulas  are  not  necessary.  Throughout  we  have  used  a method 
of  magnitude  averaging  proposed  by  Ringdal  (1976)  in  which  the  magnitudes  at 
the  individual  stations  are  assumed  to  follow  a Guassian  distribution.  Among 
this  distribution,  amplitudes  corresponding  to  some  magnitudes  will  fall  below 
the  noise  level.  Ringdal's  method  then  substitutes  a noise  measurement  at 
those  stations  which  don't  detect  and  computes  the  maximum  likelihood  estimate 
of  magnitude  based  on  measured  signals  and  noise  estimates.  The  effect  of 
this  procedure  is  to  more  accurately  define  the  magnitude  of  events  not  widely 
recorded.  The  M -m  plot  using  this  method  for  all  the  events  of  this  study 

S D 

is  shown  in  Figure  16. 


Ringdal,  F.  (1976).  Maximum-likelihood  estimate  of  event  magnitude.  Bull.  Seism. 
Soc.  Am. , 66,  789. 
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TABLE  V 


Magnitude  Data  for  East  Kazakh  Explosions 
Event  1 

Station  A M 


AVERAGE 


Event  2 


Station 


CHG 

EIL 

KON 

NAO 

AVERAGE 


Event  3 


Station 


CHG 

EIL 

KON 

NAO 

AVERACE 


Event  4 


Station 


AVERAGE 


Event  5 


Station 


NAO 

OGD 

AVERAGE 


Event  6 


Station 


MAIO 

MAO 

OGD 

AVERAGE 


*Bungum  and  TJostheim,  1976 


The  magnitude  averages  of  those  events  where  many  of  the  readings  are 
noise  levels  are  lower  than  what  would  result  if  only  measured  signal  ampli- 
tudes were  used.  Where  only  noise  measurements  were  used  in  the  M determin- 

s 

at  ion,  an  arrow  in  the  figure  represents  the  maximum  M value.  The  line  M, 

s s 

■ m^-1.5,  suggested  as  a suitable  discrimination  line  by  Blandford  and  Clark 

(1975)  and  earlier  workers,  clearly  separates  the  Baikal  events  from  the  East 

Kazakh  explosions  and  MILROW.  Note  that  Mg  for  event  6 in  our  group  is  based 

on  only  one  actual  M measurement  and  would  otherwise  lie  close  to  the  line 

s 

Vv1  .5  if  routine  magnitude  estimation  were  used.  This  line  would  also 
separate  the  suite  of  Asian  events  studied  by  Dahlman  et  al.  (1974).  The  East 
Kazakh  explosions  are  located  approximately  15  to  30  degrees  southwest  of  the 
Baikal  events. 

Corner  Frequency 

Figure  17  is  a plot  of  |Wo|  versus  corner  frequency  for  the  Baikal  events 
and  the  East  Kazakh  and  MILROW  explosions  recorded  at  LASA  and  NORSAR.  No 
long-period  P phases  were  recorded  at  the  long-period  arrays  due  to  the  small 
magnitudes  of  the  events,  so  the  values  of  |s7o|  were  estimated  from  the  spec- 
tral data  of  the  short-period  vertical  component.  Hanks  and  Thatcher  (1972) 
have  shown  that  for  a given  long-period  level  Aleutian  explosions  have  a high- 
er corner  frequency  than  Aleutian  earthquakes.  The  physical  basis  of  this 
discriminant  is  the  smaller  source-time  dimension  of  the  explosion  for  a given 
long-period  level.  However,  Figure  17  shows  no  separation  between  the  Baikal 
events  and  the  East  Kazakh  and  MILROW  explosions.  The  lack  of  a clear  discri- 
minant could  be  due  to  characteristics  of  the  sources  examined  in  this  study. 
Another  possible  explanation  is  that  due  to  the  low  magnitudes  of  the  events 


Blandford,  R.  R. , and  D.  Clark  (1975).  Variability  of  seismic  wave  forms  record- 
ed at  LASA  from  small  subregions  of  Kamchatka,  SDAC-TR-75-12,  Teledyne 
Geotech,  Alexandria,  Virginia. 

Dahlman,  0.,  H.  Israelson,  A.  Austegard,  and  G.  Hornstrom  (1974).  Definition 
and  identification  of  seismic  events  in  the  U.S.S.R.,  Bull.  Seism.  Soc. 
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Figure  17 


Long-period  spectral  level  vs.  corner  frequency  for  Baikal  events 
and  East  Kazakh  and  MILROW  explosions  from  LASA  and  NORSAR  P re- 
cordings. 


studied  here,  a clear  asymptotic  value  of  |Dq|  was  not  possible  in  most  cases. 
Long-Period  Body-Wave  Excitation 

There  were  too  few  long-period  P observations  to  determine  average  event 
ratios  of  these  quantities  to  long-period  LR  ground  displacement.  All  of  the 
long-period  P observations  were  close  to  the  noise  level,  as  would  be  expected 
for  these  small-magnitude  events.  No  clear  long-period  S waves  were  observed 
on  the  recordings  analyzed  for  the  Baikal  events. 

Depth  of  Focus 

All  of  the  Baikal  events  studied  here,  except  events  4 and  6,  had  appar- 
ent pP  phases  recorded  at  either  LASA  or  NORSAR  or  both  arrays  (see  Figure  2). 
The  pP  phase  was  not  seen  at  LASA  or  NORSAR  for  event  4;  however,  pP  for  event 
4 was  observed  at  a few  of  the  WWSSN  stations.  The  pP  phases  recorded  at  the 
WWSSN  stations  agreed  with  the  pP  depths  determined  at  LaSA  and  NORSAR  for  all 
cases.  All  the  events  were  found  to  have  crustal  foci  (Table  II).  Thus  many 
low  M^-m^  events  in  this  region  could  possibly  be  identified  as  earthquakes  on 
the  basis  of  pP  observations.  The  pP  phase  was  of  course  not  identified  for 
the  explosions,  but  lack  of  it  cannot  function  as  a positive  discriminant. 

Also  for  crustal  depth  events  there  is  very  little  pP  moveout  for  teleseismic 
stations,  so  that  pP  would  not  be  useful  as  a counter-shot-array  discriminant. 

Complexities 

We  have  computed  the  "complexity"  parameter  as  seen  at  LASA  and  NORSAR 
for  the  Baikal  events  and  East  Kazakh  and  MILROW  explosions  in  the  manner 
given  by  Lambert  et  al.  (1969).  Figure  18  shows  the  complexity  values  versus 
m,  . We  reiterate  that  LASA  and  NORSAR  beams  show  apparent  pP  signals  in  some 
cases  since  this  will  tend  to  enhance  complexity.  The  complexity  numbers  for 
the  East  Kazakh  explosions  were  low  (3.27-7.35),  as  would  be  expected  for  ex- 
plosions. The  complexity  numbers  for  the  Baikal  events  were  higher  (4.14- 
13.36),  but  overlap  the  numbers  for  the  explosions.  The  lowest  complexity 
number  for  the  Baikal  events  was  for  event  6,  a possible  explosion.  The  Bai- 
kal events  with  the  exception  of  event  6 separate  from  the  East  Kazakh  explo- 


Lambert,  D. , D.  von  Seggern,  S.  Alexander,  and  G.  Galat  (1969).  The  LONGSHOT 
Experiment,  Volume  II.  Comprehensive  Analysis,  SDL  Report  No.  234,  Tele- 
dyne Geotech,  Alexandria,  Virginia. 
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sions.  The  high  complexity  number  for  the  MILROW  signal  recorded  at  LASA  is 
possibly  due  to  a plate  refraction  effect;  LASA  is  within  the  shadow  zone  de- 
fined by  the  ray-tracing  study  of  Davies  and  Julian  (1972). 

Spectral  Ratios 

Spectral  ratios  have  been  calculated  at  LASA  and  NORSAR  according  to  a 
form  suggested  by  Lacoss  (1969) 

R = /^*55  A(f  )df  / r**  A(f  )df 

where  sums  of  equivalent  terms  of  the  discrete  Fourier  transform  have  replac- 
ed the  amplitude  spectrum  integrals  A(f)df.  These  ratios  are  plotted  against 
m^  for  the  Baikal  events  and  East  Kazakh  and  MILROW  explosions  in  Figure  19  for 

using  the  spectra  uncorrected  for  attenuation.  The  spectral  ratios  for  the 

-2 

same  P recordings  are  plotted  in  Figure  20  with  t*  values  for  a w source 

_3 

model  for  both  earthquakes  and  explosions  and  for  a w source  model  for 
earthquakes.  Bars  connect  the  spectral  ratios  for  the  two  different  earth- 
quake source  models. 

The  amplitudes  of  the  P spectra  were  greater  than  the  amplitudes  of  the 
noise  spectra  in  the  frequency  range  0.4  to  2.0  Hz  for  all  of  the  LASA  and 
NORSAR  beams.  The  spectral  ratio  does  not  apparently  depend  on  m^,  Mg,  or 
depth.  In  general  the  explosions  show  larger  spectral  ratios  than  the  Baikal 
events,  but  the  explosion  and  earthquake  populations  cannot  be  said  to  separ- 
ate. Note  that  the  spectral  ratio  value  for  event  6,  the  presumed  explosion, 
falls  along  the  border  of  the  explosion  population. 

Radiation  Pattern 

Body  waves  for  all  the  Baikal  events  except  event  6 show  both  distinct 
compressional  and  dilatational  first  motions,  suggesting  that  these  events  are 
earthquakes.  Due  to  the  low  magnitudes  of  these  events,  first  motions  were 
difficult  to  determine  on  many  records  though.  The  overlapping  compressional  * 

and  dilatational  readings  on  the  first  motion  diagrams  indicate  that  this  is 



Davies,  D.,  and  B.  R.  Julian  (1972).  A study  of  short-period  P-waves  signals 
from  Long  Shot,  Geophys.  .J . , 29,  195-202. 
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nical Note  1969-24,  Lincoln  Laboratory,  Lexington,  Massachusetts. 
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Figure  19 


Short-period  P spectral  ratio  vs.  for  Baikal  events  and  East 
Kazakh  and  MILROW  explosions  recorded  at  LASA  and  NORSAR  for  t* 
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Figure  20.  Short-period  P spectral  ratio  vs.  m.  for  Baikal  events  and  East 

Kazakh  and  MILROW  explosions  recorded  at  LASA  and  NORSAR.  The  t* 
values  are  for  w“2  source  models  for  earthquakes  and  explosions 
and  cj-3  source  models  for  earthquakes. 
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not  u reliable  discriminant.  We  were  not  able  to  determine  a radiation  pat- 
tern from  the  LR  amplitudes  owing  to  the  large  scatter,  which  is  caused  by 
proagation  effects,  such  as  Q differences,  dispersion  effects  on  time-domain 
amplitudes,  and  focussing  or  defocussing  due  to  refraction  as  illustrated  in 
the  predicted  raypaths  in  Figures  14  and  15. 


Short-Period  Excitation 


As  might  be  expected  for  these  small  magnitude  events,  there  were  no 
short-period  S observations  from  the  Baikal  events  to  determine  the  ratio  of 
S/P  ground  displacement. 


SUMMARY 


! 
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Six  events  from  the  Baikal  rift  zone,  six  East  Kazakh  explosions,  and  the 
Aleutian  explosion  MILROW  were  examined  in  a seismic  discrimination  context. 
Most  of  the  Baikal  events  lie  close  to  major  faults  observed  in  satellite  pho- 
tographs. It  was  impossible  on  the  basis  of  first  motions  to  determine  indi- 
vidual fault-plane  solutions  due  to  the  small  magnitudes  of  the  events,  but 
normal  faulting  may  be  indicated  by  the  composite  pattern.  We  were  also  unable 
to  determine  a radiation  pattern  from  the  LK  amplitudes,  probably  due  to  wide- 
ly varying  propagation  effects.  The  following  characteristics  of  the  Baikal 
events  indicate  that  all  the  events  except  event  6 are  earthquakes: 

— All  the  earthquakes  (except  # b)  had  apparent  pP  phases  recroded  at  either 
LASA,  NORSAR,  and  some  the  the  WWSSN  stations.  Although  little  moveout 
can  be  detected  for  shallow  earthquakes,  identification  of  this  phase 

is  probably  one  of  the  best  discriminants. 

— All  the  Baikal  events  studied  here  fall  above  the  line  M ” m -1.5, 

8 b 

except  event  6.  The  explosions  clearly  lie  below  the  line  M - m^-1.5. 

— Complexity  separated  Kazakh  explosions  from  Baikal  earthquakes  in  this 
data  set,  but  only  marginally  so. 

There  are  other  areas  in  Asia,  notably  near  30°N,  95°K  (Tat ham  et  al., 
1976),  where  earthquakes  do  not  clearly  separate  from  explosions  in  Asia  on 
an  plot.  Also  shot  arrays  could  be  used  to  raise  the  M^  of  an  explo- 

sion by  about  0.3  Mg  units  or  more,  placing  it  close  to  the  earthquake  pop- 
ulation studied  here. 

We  have  concluded  that  event  6 was  an  explosion  for  the  following  reasons: 

— It  was  located  north  of  the  Baikal  rift  zone  in  an  area  where  there 
are  no  active  faults.  (This  "discriminant"  would  presumably  not  be 
available  in  case  of  a test  ban  treaty.) 

— There  were  no  epicenters  from  1971  through  1975  on  the  NE1S  list  with- 
in a few  degrees  of  event  6. 


Tathara,  R.  H.,  D.  W.  Forsyth,  and  L.  R.  Sykes  (1976).  The  occurrence  of  anom- 
alous seismic  events  in  eastern  Tibet,  Geophys.  .J . , v.  45,  451-481. 
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— The  long-period  spectral  level  for  event  6 falls  within  the  explosion 
population. 

— There  are  no  apparent  pP  phases  recorded  for  event  6. 

— The  complexity  value  for  event  6 falls  within  the  Kazakh  explosion  popu- 
lation. All  the  other  Baikal  events  were  separated  from  the  Kazakh  ex- 
plosion population. 

— The  spectral  ratio  for  event  6 falls  within  the  explosion  population. 

The  M for  event  6 is  determined,  on  the  basis  of  noise  measurements  and 
8 

one  observation,  to  be  sufficiently  low  to  characterize  it  as  an  explo- 
sion. 

— An  origin  time  near  the  exact  hour  hints  at  a planned  explosion. 

In  Figure  21  we  present  SRO  recordings  of  P waves  from  this  event.  In 
each  case,  there  seems  to  be  a dilatational  first  motion.  This  contradicts 
our  conclusion  for  the  character  of  this  source,  and  we  suggest  that  these 
seismograms  are  blatant  examples  of  the  unreliability  of  first  motion  as  a 
means  of  discriminating  earthquakes  from  explosions.  The  noise  in  each  case 
might  distort  the  first  motion  considerably,  such  that  the  underlying  com- 
pressional  first  break  is  obscured. 

The  low  magnitudes  of  the  earthquakes  studied  here  (m^  from  4.48  to  5.30 

and  M from  3.56  to  5.49)  made  clear  discrimination  difficult  with  the  avail- 
s 

able  data  and  probably  represent  nearly  the  lowest  threshold  of  multi-discrim- 
inant analysis  for  the  stations  used  in  this  report.  The  installation  of  high 
quality  SRO  stations  in  Asia  will  permit  the  examination  of  lower  magnitude 
events  with  multistation  short-period  data. 

This  study  revealed  little,  if  any,  evidence  for  a highly  attenuating 
upper  mantle  zone  beneath  the  Baikal  rift.  A minor  scan  of  ISC  data  resulted 
in  no  indication  of  a generally  positive  travel-time  residual  at  stations  near 
the  zone,  in  fact  the  contrary.  The  values  of  t*  computed  from  Baikal  events 
to  NORSAR  and  LASA  were  not  significantly  larger  than  for  previous  regions 
studied  in  this  series  of  reports.  The  anomalous  zone,  if  it  exists,  must  be 
of  limited  extent,  perhaps  to  be  outlined  by  future  detailed  studies. 
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